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Abstract 

Following an electrical stimulus, the transmembrane voltage of cardiac tissue rises 
rapidly and remains at a constant value before returning to the resting value, a phe- 
nomenon known as an action potential. When the pacing rate of a periodic train of 
stimuli is increased above a critical value, the action potential undergoes a period- 
doubling bifurcation, where the resulting alternation of the action potential duration 
is known as alternans in the medical literature. Existing cardiac models treat alter- 
nans either as a smooth or as a border-collision bifurcation. However, recent experi- 
ments in paced cardiac tissue reveal that the bifurcation to alternans exhibits hybrid 
smooth/nonsmooth behaviors, which can be qualitatively described by a model of so- 
called unfolded border-collision bifurcation. In this paper, we obtain analytical solu- 
tions of the unfolded border-collision model and use it to explore the crossover between 
smooth and nonsmooth behaviors. Our analysis shows that the hybrid smooth/nonsmooth 
behavior is due to large variations in the system's properties over a small interval of 
the bifurcation parameter, providing guidance for the development of future models. 

1 Introduction 

1.1 Background 

Cardiovascular disease is the number one cause of death in the United States pQ. Over half of 
the mortality is due to sudden cardiac arrest that is often initiated by ventricular fibrillation, 
a fatal heart rhythm disorder. The induction and maintenance of ventricular fibrillation 
has been connected to the dynamics of local cardiac electrical properties [21 [3]. Therefore, 
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studying cardiac dynamics is important for understanding life-threatening arrhythmias and 
developing therapies for preventing sudden cardiac death. 

To develop an understanding of cardiac rhythm instability, we briefly review the elec- 
trophysiology of the heart. Cardiac cells respond to an electrical stimulus by eliciting an 
action potential pE], which consists of a rapid depolarization of the transmembrane voltage 
followed by a much slower repolarization process before returning to the resting value (Fig. 
[TJ. The time interval during which the voltage is elevated is called the action potential 
duration (APD). As shown in Fig. [TJ the time between the end of an action potential to the 
beginning of the next one is called the diastolic interval (DI). The time interval between two 
consecutive stimuli is called the basic cycle length (BCL). 
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Figure 1: Schematic action potential showing the response of the transmembrane voltage to 
periodic electrical stimuli. 

Under a periodic train of electrical stimuli, the steady-state response consists of phase- 
locked action potentials, where each stimulus gives rise to an identical action potential (1:1 
pattern) when the pacing rate is slow. When the pacing rate becomes sufficiently fast, 
the 1:1 pattern may be replaced by a 2:2 pattern, so-called electrical alternans [51 EE E], 
where the APD alternates between short and long values. Using theory and experiments, a 
causal connection between alternans and the vulnerability to fatal cardiac arrhythmias such 
as ventricular fibrillation has been established by various authors [2, 3, 5-13]. Therefore, 
understanding mechanism of alternans is a crucial step in detection and prevention of fatal 
arrhythmias. 

It has long been hypothesized [HI El QUI Ell E21 E3j that alternans is mediated by a classical 
period-doubling bifurcation, which can be described using a smooth iterated map, and which 
occurs when one eigenvalue of the Jacobian crosses the unit circle through —1 [TJ1 E5]- We 
restrict our attention to supercritical rather than subcritical bifurcations because the former 
are observed in most experiments and theoretical models exhibiting electrical alternans. 
Based on this hypothesis, various authors attempted to develop criteria for the onset of 
alternans [HI EH EU] as well as algorithms to control alternans [UJ E2J E_3] . Recently, a few 
authors [161 E3 EH1 EH1 ED] proposed a different hypothesis: alternans may be mediated 
through a border-collision period-doubling bifurcation. Border-collision bifurcations occur 
in piecewise smooth maps [2TI [22] . In contrast to classical period-doubling bifurcations, 
eigenvalues are not indicative of the onset of a border-collision period-doubling bifurcation. 
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Instead, a border-collision bifurcation occurs when a branch of fixed points collides with 
a border, i.e., a discontinuity surface in state space. Knowing the mechanism of alternans 
may help researchers to choose the proper types of functions to model this instability. More 
importantly, to develop model-based control methods requires knowledge of the underlying 
dynamics [T7j. 

The aforementioned intrinsic differences between the two bifurcation types lead to dif- 
ferences in their bifurcation diagrams, as depicted in Fig. [2j Here, the two bifurcated 
branches of a smooth period-doubling bifurcation become tangent to each other at the bifur- 
cation point, while the bifurcated branches of a border-collision bifurcation open at an angle. 
Thus, in principle the bifurcation diagrams should distinguish between the two bifurcation 
types. However, in practice, experiments can provide only a limited number of measure- 
ments (especially in biological systems), so the resulting bifurcation diagrams do not have 
sufficient resolution. This is illustrated in Fig. [21 where the discrete points representing 
experimental data along a bifurcation diagram do not readily reveal the true type of bifur- 
cation. Therefore, there is a need for a more sensitive technique to differentiate between the 
two bifurcations. 

(b) 

>- 

i 

B hi{ B B hi{ B 

Figure 2: Schematic bifurcation diagrams of period-doubling bifurcation: (a) a smooth type 
and (b) a border-collision type. Here, B represents a bifurcation parameter and A represents 
fixed-point solutions. 




1.2 Prebifurcation Amplification 

Based on prebifurcation amplification, our group has developed a robust technique to distin- 
guish between smooth and border-collision bifurcations [27-30]. Here, we briefly review the 
results. It has been shown theoretically and experimentally that, near the onset of a smooth 
period-doubling bifurcation, subharmonic perturbations in a bifurcation parameter result in 
amplified disturbances in the response, a phenomenon known as prebifurcation amplification 
[23| [2H 123 [2H]. In the following, we will show that, under variations in system parameters, 
prebifurcation amplification exhibits qualitatively different scaling laws in border-collision 
period-doubling and smooth period-doubling bifurcations. Thus, prebifurcation amplifica- 
tion is a useful technique to distinguish between the two possible types of period-doubling 
bifurcations. 

To illustrate the concept of prebifurcation amplification, we consider a dynamical system 
described by the following map 

x n+1 = f (x n ; B) , (1) 
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where B represents a bifurcation parameter, e.g., the BCL in cardiac models. Both the 
function / and state variable x may be one- or multi-dimensional. Let us assume that, at a 
critical value B = Buf, the system undergoes a period-doubling bifurcation that is either a 
smooth type for smooth / [T4l 115] or a border-collision type for piecewise smooth / [211 122] . 
We further assume that the stable period-one solution lies on the side B > Buu as indicated 
in Fig. [3 

When a subharmonic perturbation is applied to B under conditions when B > B^, it 
renders map ([I]) as 

x n+1 = f(x n ;B + (-l) n 5), (2) 

where 5 is the amplitude of the perturbation. The perturbation may also be imposed in the 
form of B — (—1)™ 5, which leads to a solution only different in phase from that of Eqn. ([2]). 
Since B represents the pacing interval in cardiac models, such a variation in B is referred 
to as alternate pacing. For B greater than but close to and small 5, the steady-state 
response of Eqn. (j2J) consists of alternating recurrent states of x even and x dd, which satisfy 
the following conditions 

aWen = / (awd; B -5) , (3) 

^odd = / Oeven; B + 5) . (4) 

In cardiac models, one component of the vector x is the APD, henceforth denoted by A. 
Alternate pacing of these models results in a long-short beat-to-beat variation in pacing 
intervals, which in turns cause alternation in A even when B > B^. Since a period-doubling 
bifurcation is sensitive to subharmonic perturbations, perturbations in B result in amplified 
disturbances in A. The effect of prebifurcation amplification can then be characterized by a 
gain defined as follows 

14 — A , J 

-p | ^even ^-odd | 

= 25 ' 
1.2.1 Gain of Smooth Bifurcations 



(5) 



Several authors [2_2[23],[25j[26] have investigated the influence of parameters on prebifurcation 
amplification in smooth period-doubling bifurcations. In a previous paper [28] , we explored 
the scaling laws between the amplification gain T and the parameters B and 5, using a 
mapping model of arbitrary dimension. It was shown there that the gain of a smooth 
bifurcation satisfies the following relation, 

c5 2 T 3 + (B-B hi{ )T- \k\ =0, (6) 

where c and k are constants determined by the system's properties at the bifurcation point. 
It was established that the gain is infinite if and only if B = B^ and 5 = 0. The rate 
of divergence as the parameters tend to (B hi{ , 0) depends on the path taken. For example, 
when 5 is extremely small, the gain tends to infinity as [B — B^)^ 1 ; on the other hand, 
when B = B^, the gain tends to infinity as 5~~ 2 / 3 . 
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In cardiac experiments, it is very difficult to accurately locate the bifurcation point. 
Moreover, the existence of noise and the limitation on the number of measurements restrict 
one from using very small perturbations. Instead, one can investigate the gains under two 
protocols: i) let B approach B^f while retaining a finite and constant 5; and ii) let 5 approach 
zero while retaining a constant B > Bu{- As has been established in [28] that, under constant 
5, T scales according to (B — B^t) 1 except when B — B^i is sufficiently small, where the gain 
becomes saturated. Alternatively, under constant B > B^, T scales to <5~ 2 / 3 except when 5 
is sufficiently small, where the gain becomes saturated. Figure [3] (a) and (b) schematically 
show the behaviors of T vs. B and T vs. 5, respectively. 



1.2.2 Gain of Border-Collision Bifurcations 

The system ([1]) possesses a border-collision bifurcation if the function / is piecewise smooth 
as follows 

f („. f>\ _ / h ( x ; B ) > if h ( x ) < (7 \ 

J[X,B) | f 2 ( X ;B), ]fh(x) >0 ' {<) 

where h is a smooth scalar function and h (x) = indicates a "border" in the state space, 
on which fi (x; B) = /2 (x;B). An approximate expression for the gain is derived in the 
Appendix using a one-dimensional map; results for general maps can be found in [29|. To 
lowest order, the gain is piecewise smooth as follows 

J- const) 

if {B - B hif ) /5 > p cr it / Q x 

Tconst - 7 (2=f* - Pcrit) , if (B - B hii ) /8 < Pcrit ' 1 ' 



where r const , 7, and p cri t are positive constants determined by system properties. Therefore, 
the gain is a constant along any straight line (B — £?bif) /8 = const. Since all these lines 
intersect at (B hii , 0), the gain at this point is not defined. 

Again, we apply the two protocols described in the previous subsection. When 5 is 
constant, the gain is constant when B > B crit = B^f + p cr it^ an d varies linearly as B 
when B < B CT i t . Alternatively, when B is constant, the gain is constant when 5 < 5 cr it = 
(B — £>bif) /pcrit and varies as <5 -1 when 5 > 5 CT i t . Schematic diagrams of these behaviors are 
shown in Figs. [3] (c) and (d). 

It is evident from Fig. [3] that behaviors of the gain are qualitatively different for a 
smooth bifurcation and for a border-collision bifurcation. However, the differences between 
Figs. [3] (a) and (c) may be difficult to detect for discrete data or for data disturbed by 
noise. Conversely, differences in Figs. [3] (b) and (d) are apparent even for discrete data 
and in the presence of noise. Therefore, investigating the T vs. 5 under alternate pacing 
provides an unambiguous way to distinguish between the two bifurcations. Moreover, since 
this technique relies on the trend of the gain rather than the magnitude, it allows one to 
distinguish between smooth and nonsmooth behaviors in experiments without the need to 
accurately locate the bifurcation point 
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Figure 3: Prebifurcation gain T of a classical period-doubling bifurcation (a,b) and of a 
border-collision period-doubling bifurcation (c,d). In panels (a) and (c), 5 stays constant; 
in panels (b) and (d), i?bif < B stays constant. Comparison between panels (b) and (d) 
provides the most revealing difference between the two bifurcations. 

1.3 Hybrid Behavior of the Prebifurcation Gain 

To identify the bifurcation mechanism mediating cardiac alternans, we implemented the 
aforementioned technique in paced in vitro bullfrog heart [30l [3T] , where the experiments 
reveal a novel phenomenon that cannot be explained by the above simple dichotomy of 
smooth/nonsmooth bifurcations. Specifically, our experiments show that very close to the 
bifurcation point, V decreases with 5, which agrees with the smooth bifurcation (Fig. 3b), 
whereas further away T increases with 5, which agrees with the border-collision bifurcation 
(Fig. 3d). A bifurcation that exibits such a crossover between smooth and border-collision 
behaviors is named a hybrid period-doubling bifurcation [301 13T] . We further found that the 
essence of this hybrid behavior can be reproduced by a model of a so-called unfolded border- 
collision bifurcation. In the remainder of this paper, we will carry out a detailed analysis of 
the unfolded border-collision model. This analysis will help to understand the mathematical 
mechanism underlying the crossover between smooth and nonsmooth behaviors, providing 
guidance for the development of future models. 

2 A Model of an Unfolded Border-Collision Bifurca- 
tion 

We explore the mechanism of the aforementioned hybrid behavior using the unfolded border- 
collision model presented in [301 EI] • F° r illustration purpose, we first consider a piecewise 
smooth map 

A n+1 = A c + a (D n - Ah) +P\D n - Ah | , (9) 

where A n and D n denote the nth action potential duration and diastolic interval, respectively. 
Note that D n = B — A n as can be seen from Fig. [TJ Under the following conditions (cf. 
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[22]) 

-l<a + /?<l<a-/?and - 1 < a 2 - /? 2 < 1, (10) 

map ([S]) possesses a border-collision period-doubling bifurcation at 

B C = A C + D th . (11) 

To remove the nonsmoothness of map Qj, we "unfold" the singular term (3 \D n — Ah I as 
follows 



A n+1 = A c + a (D n - Ah) + P\f (A. - Ah) 2 + A 2 - (12) 

Map (fT2l) represents a one-parameter family of maps that reduces to map (J9]) when A = 0. 
For any A 7^ 0, the unfolded map (Tl2|) is smooth and exhibits what is technically a smooth 
period-doubling bifurcation. Nevertheless, the dynamics of map (Q and map (fl2j) exhibit 
no significant differences except when B — B c is less than or on the order of A- It is worth 
noting that there are other ways to unfold the border-collision map (JH]). Here, we choose 
map (fT2l) because of its simplicity and ease of analysis. 

In the following, we show that map ffl2|) has a smooth period-doubling bifurcation if 
A 7^ 0. To this end, we denote the bifurcation point by A = A hi f = Aif — Aif and let the 
Jacobian of map ( fl2l) equal to —1 at the bifurcation point; in symbols 

g (Aif - Ah) . 

— a = = —1. (13) 

y (Aif - Ah) 2 + A 2 

It follows from Eqn. (fl~3l) that 

/? (Aif - Ah) = (1 - «) \/ (Aif- Ah) 2 + A 2 - (14) 

Since /3 < as can be shown from the conditions in Eqn. ffTU]) . the term Dbif — Ah has an 
opposite sign as the term 1 — a; in symbols 

(Aif- Ah)(i-«) <o. (15) 

Evaluating Aif from Eqn. f[T4l) and considering the conditions ffTUl) and ffTBT) yields 

(l-a)A 

Aif = Ah - / ; = ■ (16) 
yJP - (1 - a) 2 

Thus, APD at the bifurcation point can be written as 



A hi{ = A c + a (Aif - Ah) + N (Aif - Ah) 2 + A 2 



S ' 



A c + a (Aif - Ah) H : , (17) 

1 — a 

jf? - (1 - af 
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and the corresponding value of BCL is 



Bbif = ^bif + -Dbifj 



A c + Ah - 



(l-a 2 + (3 2 )D s 
sJJ 2 " (i - «) 2 



(18) 



Comparing I?bif and B c reveals that the smooth period-doubling bifurcation in map (fT2l 
reduces to the border-collision period-doubling bifurcation in map as D s —>■ 0. Moreover, 
it can be shown from Eqn. (fTUj) that 1 — a 2 + (3 2 > so that < B c . Figure H] demon- 
strates schematically the relation between a border-collision bifurcation and the unfolded 
bifurcation. 



Figure 4: Schematic diagram showing a border-collision bifurcation (solid) and the unfolded 
bifurcation (dashed). 



2.1 Analysis of the Response to Alternate Pacing 

To study the prebifurcation amplification of map (Tl2|) . we apply an alternating perturbation 
to the BCL's; in symbol, B n — B + (—1)™ 6, where B is a baseline BCL and 5 is a small but 
nonzero perturbation. Under this alternate pacing, it follows that D n — B + (— l) n 5 — A n . 
Here, we require that B > B^ because prebifurcation dynamics is of interest. Denoting the 
steady-state APDs under alternate pacing by A even and A odd , it follows from Eqn. (Tl2|) that 



A 




B 




(19) 



(20) 



where 



even 



B + 5-A, 



B-8-A, 



-even • 



odd; 



(21) 
(22) 
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For later convenience, we let 

B = B c + AB = A c + Ah + AB (23) 
and we define A even and A odd by 

A even = A even + A c and A odd = A odd + A c . (24) 
Substituting the above equations into Eqns. (fT9l and (1201) yields 



A even + a (A odd + 5 - AB) = (3\/(A odd + 5 - AB) 2 + B s 2 , (25) 



A odd + a (A evcn - 5 - AB) = (3y/(A eYen - 5 - AB) 2 + D 2 S . (26) 

One can then show that 

((1 - a) (A cvcn - A odd ) + 2 a 5) ((1 + a) (A cvcn - A odd ) -2aAB) 

= f3 2 (A odd - A even + 2 5) (A odd + A cvcn - 2 AB) (27) 

Let 

A even A odd /r>o\ 

7 = 25 ; (28) 

i.e., 7 is a gain-like quantity that can be either positive or negative (cf. [5]). Substituting the 
definition of 7 into Eqn. (|27|) yields 

2 5 ((1 - a) 7 + a) ((1 + a) (A cvcn - A odd ) - 2 a AB) 

= 25 (3 2 (1 - 7) (A odd + A cvcn — 2 AB) . (29) 

Because we consider nonzero 5, Eqn. (I2"9l can be reduced to 

7 ( Cl AB + c?! (A even + A odd )) = c 2 AB + d 2 (A even + A odd ) , (30) 

where 

ci = -2(/3 2 + a(l-a)), (31) 

dj = 1 - a 2 + /5 2 , (32) 

c 2 = 2 (« 2 - (3 2 ) , (33) 

d 2 = - (a (1 + a) - f3 2 ) . (34) 



It then follows that 



= c 2 AB + rf 2 (A even + A odd ) 
c\AB + d a (A cvcn + A odd ) ' 



Since A ev en and A odd depend on AB and 5, 7 is a function of AB and 5. 
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Recalling definition (J5J), we find the prebifurcation amplification gain as 



c 2 AB + d 2 (A 

even + A odd ) 



c x AB + di (A even + A odd ) 
Particularly, when AB = 0, i.e. B = B c = A c + D^, the gain is 



d 2 



a + (a 2 - P 2 ) 
1 - {a 2 - (3 2 ) 



(36) 



(37) 



Therefore, when B = B c , V is same for all 5. With some manipulation, one can show that 
dT/dB at B = B c . Moreover, when B is sufficiently close to Buf < B c and 8 is fixed, 
T decreases as B increases as described in previous section and proven in [28]. Thus, for a 
given 5, T is a monotonically decreasing function of B and T becomes constant at B = B c . 
Because T is a monotonically decreasing function of 5 when B > -Bbif, as shown in the 
previous section (see also [28J), it follows by continuity that T will increase as 5 increases in 
the region of B > B c . In other words, the map ffl2|) exhibits smooth like behavior when B 
is sufficiently close to -Bbif and border-collision like behavior when B > B c (see the relation 
between Buf and B c in Fig. H]). 



2.2 Numerical Example 

Before comparing the proposed model to experimental data, we review the class of models 
that are most commonly used in the cardiac research community. These models relate APD 
and DI through exponential functions. Typically, parameters of a model are obtained by 
fitting the model to the so-called dynamic restitution curve, which is a plot of the steady- 
state APD vs. DI. For example, in their pioneering work, Guevara et al. [32J proposed a 
model of cardiac dynamics as 

A n+1 = 201 - 98 e - Dn/i3 - 35 e - Dn/653 , (38) 

where all variables and parameters have the unit of millisecond. Parameters of map (|38l) 
were obtained by fitting the dynamic restitution curve measured in experiments performed on 
quiescent aggregates of ventricular cells from 7-day-old embryonic chick hearts [32] • Although 
the model of Guevara et al. fits the dynamic restitution curve reasonably well (Fig. \5\ top 
panel), it does not accurately describe the response beyond the bifurcation to alternans, as is 
evident from the bifurcation diagram of steady-state APD vs. BCL (Fig. [5l bottom panel). 
A careful examination reveals that the dynamic restitution curve is well approximated by two 
distinct parts with significantly different slopes. The transition between the two slopes occurs 
with a small interval near DI^ 60 ms, which is also approximately where the transition to 
alternans occurs. Now, we recall that the unfolded border-collision model (JEJ), with properly 
chosen parameters, describes such rapid changes between two distinct slopes. Fitting map 
( TT2|) to the experimental dynamic restitution data, we obtain a set of parameters 

a = 0.69, (3 = -0.64, 
A c = 161ms, D t h = 62 ms, and D s = 15 ms. (39) 
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As shown in Fig. [5j the unfolded border-collision map (|T2|) with these parameters faithfully 
reproduces the bifurcation diagram, including the alternans branches. As demonstrated in 
Fig. HI the bifurcation diagram of a border-collision map is close to that of its unfolded 
counterpart except near the bifurcation point. Thus, one expects a reasonable fit to the 
experimental data in Fig. [5] using a pure border-collision map, i.e., letting D s = 0. How- 
ever, as has been established in the previous section, the nonsmooth map can not capture 
hybrid behaviors in the prebifurcation gain. Here, for clarity, the bifurcation diagram of the 
corresponding border-collision model is not shown in Fig. [5j 
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Figure 5: Comparison between the model of Guevara et al. (solid) and the unfolded border- 
collision model (dashed) in fitting the experimental data in [32] (points). Although both 
models fit the dynamic restitution curve well (top panel), the unfolded border-collision model 
fits alternans data much better (bottom panel). 

We then simulate map (|T2j) with alternate pacing. Figure [6] shows V vs. B for different 
values of 5. These curves cross one another at B = B c = 223 ms. Note that a period-doubling 
bifurcation occurs at B = 198 ms. It is clear that T vs. 5 displays a trend consistent with 
a smooth bifurcation (cf. Fig. [3] (b)) when B < B c and, on the other hand, T vs. 5 
shows a trend consistent with a border-collision bifurcation when B > B c (cf. Fig. [3] (d)). 
Since Guevara et al. did not perform alternate pacing experiments, no data are available for 
comparison. However, we note that the simulation here is in qualitative agreement with our 
previous experiments on bullfrog ventricles [301 EI]- 
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Figure 6: Prebifurcation amplification predicted by the unfolded border-collision model ( |T2l) : 
5 = 5 ms (solid), 5 = 10 ms (dashed), and 5 = 15 ms (dotted). 

3 Discussion and Conclusion 

Theoretical analysis of the prebifurcation amplification reveals that different scaling laws are 
associated with smooth and border-collision period-doubling bifurcations. The differences 
appear in the following three aspects. First, the gain of a smooth bifurcation tends to infinity 
as (B,S) approaches (-Bbif, 0); conversely, the gain of a border-collision bifurcation is finite 
everywhere but not defined at (-Bbif, 0). Second, the gain of a smooth bifurcation varies 
smoothly under changes in system parameters while that of a border-collision bifurcation 
undergoes a nonsmooth variation as parameters cross a boundary in the parameter space 
(see Figs. |3] (a,b) and [3] (c,d)). Third, under constant B and increasing 5, the gain of a 
smooth bifurcation decreases while that of a border-collision bifurcation increases. Thus, 
the gain versus perturbation size relation provides a more sensitive criterion to differentiate 
between the two bifurcation types. As can be seen from Fig. [3] (b) and (d), even with few 
data points, the T vs. 5 relation clearly reveals the underlying bifurcation mechanism. On 
the other hand, the bifurcation diagram does not allow one to distinguish between the two 
bifurcations with only a few data points nor does the T vs. B relation. 

Although the technique described here was developed with a goal of identifying the type 
of bifurcation mediating alternans, the analysis is based on general iterated maps. Thus, 
the results are independent of any physical details of cardiac dynamics and can be readily 
applied to any dynamical systems. 

The analysis based on simple dichotomy of smooth/border-collision bifurcation has lim- 
itations. Since it is assumed that a system either has well behaved derivatives or is discon- 
tinuous in first derivatives, the result is not directly applicable to the intermediate case, i.e., 
a system whose first derivatives are continuous but change rapidly. The model of unfolded 
border-collision bifurcation studied here serves to address the latter case. 

Previous experimental findings [30, 31J suggest that modeling of cardiac dynamics should 
consider the rapid changes in the system's properties, i.e., large variations over a narrow 
parameter interval. As one example, we study here the smoothed version of a border- 
collision model. We show that the smoothed map indeed unfolds the original border-collision 
period-doubling bifurcation to a smooth one. In addition, we carry out the analysis of the 
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unfolded map under alternate pacing. The result indicates that the unfolded border-collision 
model exhibits hybrid smooth/nonsmooth behaviors, which is in qualitative agreement with 
previous experimental observations on bullfrog hearts [30j [31] . We further illustrate that 
the unfolded border-collision model can more accurately describe alternans observed in an 
experiment on embryonic chick hearts [32]. The fact that hybrid behaviors are observed 
in different species indicates that this phenomenon may be prevalent in cardiac dynamics. 
It is worth noting that, besides the model studied here, the crossover between smooth and 
nonsmooth behaviors can also be captured by other types of maps. We choose the current 
model solely based on its simplicity and ease of analysis. 

We note that many other physical systems also possess rapid changes in systems' prop- 
erties. To fully describe such rapid changes, one would need to use functions with highly 
localized properties. For convenience of analysis, these highly localized functions are often 
replaced by piecewise smooth functions, where each piece adopts a much simpler form. Per- 
haps, the simplest example is a bouncing ball, whose velocity changes rapidly before and 
after impacts and is often modeled by an instantaneous jump using the coefficient of resti- 
tution (see other examples in a recent special issue of the journal Nonlinear Dynamics on 
discontinuous dynamical systems [33J). Although this approach has proven to be useful in 
many problems in engineering and science, it brings up a more subtle question on the rela- 
tion between the piecewise smooth bifurcation problem and the original smooth bifurcation 
problem. In [M], Dankowicz purposefully coarsened a smooth vector field with a piecewise 
smooth one and compared their bifurcation diagrams. The full potential and limitations of 
the idea of intentional nonsmoothing of a smooth function need to be explored in future 
research. 
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Appendix: Alternate Pacing of a Border-Collision Map 

In a previous paper [29], we have shown the general results of prebifurcation amplification 
for border-collision bifurcations using high-dimensional maps. Here, we briefly review the 
results using a one-dimensional map for simplicity. Consider a one-dimensional piecewise 
continuous map of A with a bifurcation parameter B as follows 



where fi (A; B) = f 2 (A; B) when A = A hii . Assume a border-collision bifurcation occurs at 
B = Sbif and A = A^, as indicated in Fig. [2](b). Then the following conditions are satisfied 
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at the bifurcation point 

d A f 2 < -1< d A h < 1, (41) 
< d B fi = 9b h = d B f, (42) 

where all derivatives are evaluated at the bifurcation point (Abif; -Bbif)- For conditions on 
border-collision period-doubling bifurcations in mult i- dimensional maps, see [2T1 129] . 
Alternate pacing changes the map (|40p to 



/ / 1 (A n ; J B + (-l) n 5), ifA n >A bi f 
n+1 I / 2 (A»;S + (-l) n 5), ifA n <A bi f ' l4dJ 

Due to the alternating perturbation, the steady state of Eqn. fj43|) is a period-two solution, 
whose two branches can be written as 

A _ f A od d (5, 5) , for odd ra ^ , 

\ A even > f° r even n 

Particularly, 

Add {Buu 0) = A eve n {B uh 0) = A bif . (45) 

This solution consists of two different types: 1) in a unilateral solution, both branches are 
above the border, i.e., A even > A hii and A odd > A bif ; 2) in a bilateral solution, one branch is 
above and the other branch below the border, i.e., (A even — Aut) X (A odd — A hi{ ) < 0. In the 
following, we restrict attention to B > B^ (prebifurcation condition) and deal with the two 
types of solutions, respectively. 

Unilateral Solution 

Because Aeven > Abif and A Q dd > Abif, it follows from Eqn. ( T43l) that 

Aeven = f\ (A odd] B - 5) } (46) 
Aodd = /l (Aeven! B + 5). (47) 

To leading order, the solution of Eqs. (JIB]) and fj47|) is 

Aeven = A bif + - 9b / (B - B m ) - sf 8, (48) 

Aodd = Abif + : (B - B m ) + - d * f 5. (49) 

1 - oaJx 1 + o A fi 

Recalling the conditions in Eqs. fl4"Tj) and (jUJ), it follows that A Q dd > A even . Moreover, it 
follows from Eqn. (148]) that the unilateral solution is valid as long as 

B - B m > p ciit 5, (50) 

where 

1 - d A fi „ 

Pcrit = — > 0. 51 
1 + o A fi 
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Bilateral Solution 

The bilateral solution occurs in the region B — < p cr - lt 5. By continuity, the solution in 
this region satisfies A odd > Auf > A evcn . It follows from Eqn. ( j4"3l) that 



A even = fx (A odd ; B - 5) , (52) 
A odd = f 2 (A cvcn ,B + 5). (53) 

Linearizing Eqns. (1521) and fl53|) around A = v4 bif and B = B hi f yields 

A evcn = A hii + d A h * (A odd - A bif ) + d B f * (B - B Ui - 5) , (54) 
^4odd = A hii + d A f 2 * (A cvcn - A hii ) + d B f * (B - B ui + 6) , (55) 

where the derivatives are evaluated at (A^a ; £>bif) • Solving the above equations yields the 
leading-order solution for A even and A odd as 

A A i + d A fi) (B - B hii ) - {\-d A h)5 a 

A ev en = ^bif H : a , Q , <W, (56) 

4 4 ■ (1 + d A f2) (B - But) + (I - d A f 2 ) 5 a 

A odd = A hi{ H — — — d B j. (57) 

1 - <W 2 Oa/i 

Prebifurcation Gain 

When B — £> bif > p crit 5, it follows from Eqs. and (I4"9"j) that the gain is 

r _ A}dd — ^even _ 9 B f _ p 

1 - 26 - ttoj; = const - (58) 

When B — £> bif < p crit 5, it follows from Eqs. (156"]) and (!57|) that the gain is 

^4odd A even /rrv^ 

1 = 2l (59) 

= r const - 7 f — — - p cvlt J , (60) 



where 
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